Systematic Study of Fission Barriers of Excited Superheavy Nuclei 
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A systematic study of fission-barrier dependence on excitation energy lias been performed using 
the self-consistent finite-temperature Hartree-Fock+BCS (FT-HIF-I-BCS) formalism with the SkM* 
Skyrme energy density functional. The calculations have been carried out for even-even superheavy 
nuclei with Z ranging between 110 and 124. For an accurate description of fission pathways, the 
effects of triaxial and refiection-asymmetric degrees of freedom have been fully incorporated. Our 
survey demonstrates that the dependence of isentropic fission barriers on excitation energy changes 
rapidly with particle number, pointing to the importance of shell effects even at large excitation 
energies characteristic of compound nuclei. The fastest decrease of fission barriers with excitation 
energy is predicted for deformed nuclei around A'^=164 and spherical nuclei around A''=184 that 
are strongly stabilized by ground-state shell effects. For nuclei ^^''Pu and ^^®Fm, which exhibit 
asymmetric spontaneous fission, our calculations predict a transition to symmetric fission at high 
excitation energies due to the thermal quenching of static reflection asymmetric deformations. 

PACS numbers: 24.75.-|-i, 21.60.Jz, 27.90.+b, 24.10. Pa 
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I. INTRODUCTION 

The mere existence of the heaviest and superheavy 
nuclei with Z>10i is primarily determined by shell ef- 
fects [H, H, S i, i, i, I3- The ground-state (g.s.) shell 
corrections also determine fission barriers of those sys- 
tems d, [m as their liquid-drop fission barri- 
ers are negligible. The discoveries of new elements using 
the cold- and hot-fusion reactions over the last 
decade provide us with fundamental information about 
the structure of the nucleus and the possible existence of 
the "island of stability" at the limit of the nuclear mass 
and charge. 

Since the cross sections for production of superheavy 
nuclei using combinations of available stable projectiles 
and targets are exceedingly low, the major experimental 
challenge is to find optimal conditions that would lead 
to the synthesis of the species of interest [l^ [ij, [l5| . 
Isotopes of elements with Z up to 113 have been pro- 
duced in cold-fusion reactions using lead or bismuth tar- 
gets. In these experiments, the compound nucleus (CN) 
is formed at relatively low excitation energies E* of '^lO- 
12 MeV. Recently, using the beams of ^^Ca and actinide 
targets, superheavy elements with Z=112-116 and 118 
have been synthesized lj| . The compound nuclei formed 
in such hot-fusion reactions are more neutron-rich than 
those produced in cold- fusion experiments, and they are 
significantly more excited, i?*~36-40 MeV. 

The crucial quantity that determines the synthesis of 
superheavy elements is the CN survival probability [Tol . 
llSl . [TgI . [TtI , which strongly depends on the fission barrier 
characteristics. Since shell effects arc quenched at high 
temperatures (see, e.g., Refs. [iS, 19, 20, 21, 22, 23]), the 
stability of the heaviest and superheavy elements with 
respect to particle emission and fission is expected to 



strongly depend on excitation energy. 

In the previous paper (2^ . it was demonstrated that 
fission barriers of excited superheavy nuclei vary rapidly 
with particle number. The main objective of the present 
study is to address this question globally by performing 
systematic calculations of fission barriers of superheavy 
nuclei as a function of excitation energy Our survey has 
been carried out within the nuclear density functional 
theory (DFT) generalized to finite temperatures. Guided 
by results of Ref. [l^], we assume that the fission pro- 
cess is isentropic in character. The effects due to the E* 
dependence of triaxial and reflection asymmetric defor- 
mations are quantified and the resulting barrier damping 
parameters are extracted. 

We also investigate the transition from asymmetric to 
symmetric fission with increasing excitation energy. Ex- 
perimental studies indicate that there is a system- 
atic increase in the symmetric mass yield relative to the 
asymmetric one with excitation energy. By calculating 
the refiection-asymmetric deformations along static fis- 
sion pathways, we show that such a transition indeed 
takes place in selected nuclei. 

The manuscript is organized as follows. Section |TT] 
briefly summarizes the FT-HFB formalism. In partic- 
ular, the need for an isentropic, rather than an isother- 
mal, description of the fission process at finite excitation 
energy is emphasized. The particular realization of the 
FT-HF-I-BCS model applied in our work is presented in 
Sec. mil Excitation-energy dependence of fission path- 
ways for two representative nuclei, ^^°Pu and ^^®Fm, is 
discussed in Sec. IIVI together with the results of our sys- 
tematic calculations of the excitation-energy dependence 
of the inner fission barrier of superheavy elements. Our 
survey clearly demonstrates that the damping of the first 
barrier with E* exhibits an appreciable dependence on 
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shell effects. Finally, the summary of our work is con- 
tained in Sec. [V] 



II. FINITE-TEMPERATURE HFB APPROACH 



fission path. In this way, the temperature changes with 
deformation. In particular, the temperature of the lowest 
minimum is always greater than that of the first barrier, 
and this difference is crucial for the fission barrier damp- 
ing. 



Within the mean-field approach, heated nuclei can be 
self-consistently treated by the finite-temperature DFT, 
either within Hartree-Fock (HF) [H, [13, [H, [ll] or, if 
pairing is considered, in the Finite- Temperature Hartree- 
Fock-Bogohubov (HFB) method [2a, .30, 31,^12, 33]. The 
equilibrium state of a nucleus at a fixed temperature T 
and chemical potential /i is obtained from the minimiza- 
tion of the grand canonical potential 0, [s^l : 



E -TS - fiN, 



(1) 



where E=Tr{DH) is 
S=-kTT{D\nD) is the 
particle-number, and the density operator D given by 



the average energy, 
entropy, N=Tt{DN) is the 



(2) 



with (3 = 1/kT. In the mean-field approximation, the 
two-body density operator defined in Eq. ([2]) is replaced 
by a one-body counterpart. The variation of Q, with 
respect to density, leads to the temperature-dependent 
HFB equations [35|: 
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(3) 



where 7i is the temperature-dependent HFB Hamilto- 
nian. Finite-temperature particle and pairing density 
matrices [3^ in the FT-HFB formalism are given by 



p{f3) = UfU^ + V*{l-f)V, 
= UfV^ + V*{1~ f)U, 



(4) 
(5) 



and depend on the Fermi occupations = (1 + e'' 

The isothermal scenario, sometimes assumed in the 
context of fission process [2^ 113], cannot be correct as 
the compound nucleus is not in contact with a heat bath. 
Considering the fission as an adiabatic process, the isen- 
tropic picture seems to be more appropriate [Ij, |3^ . As 
discussed in Refs. [l^, [13, [11], the two descriptions of 
fission can be operationally related through the thermo- 

which simply 



dynamical identity ^- 
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\dQ2O )r 
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states that the generalized driving force associated with 
the deformation Q20 depends only on the state of the sys- 
tem. This identity, useful in practical calculations, has 
recently been verified numerically in Ref. ^24] wherein the 
importance of self-consistency has been pointed out. 

In this work, we shall follow the isentropic picture. The 
entropy S = S{T) has been defined as in [2^], i.e., it 
corresponds to the free energy minimum at temperature 



III. THE MODEL 

Barrier heights obtained within the HFB and HF+BCS 
approaches are quite similar at low temperatures [37.. JiSl] . 
Moreover, beyond kT ~ 0.7 MeV, the two approaches 
are identical as the static pairing vanishes [23, [H, [33 ]. 
For that reason, in this study we shall present the FT- 
HF-hBCS results only 

Our FT-HF-I-BCS calculations were carried out with 
the Skyrme SkM* functional [1^ in the particle-hole 
channel. This functional has been optimized at large 
deformations; hence, it is often used for fission barrier 
predictions. In the pairing channel, we employed the 
density-dependent delta interaction in the mixed variant 

m-. 



Vir - r') = ^0 (1 - p(r)/2po) <5(r - r') 



(6) 



T. This value of S is then kept fixed along the 



where po = O.lGfm^^. The pairing-active space in BCS 
was assumed to consist of the lowest Z/N proton/neutron 
HF levels. The pairing interaction strengths Vq are —438 
and —372 (in MeVfm'^) for protons and neutrons, respec- 
tively. They were adjusted to reproduce the experimental 
odd-even mass differences in ^^^Fm. 

It is known from numerous studies [^, [ll|, [H, [HI that 
the first saddle point is lowered by several MeV by tri- 
axial degrees of freedom and that beyond the first bar- 
rier refiection-asymmetric deformations may become im- 
portant. Therefore, when studying saddle points and 
fission pathways, it is imperative to employ a model 
which is capable of breaking axial and mirror symme- 
tries simultaneously. For that reason, we employed a 
symmetry-unrestricted DFT solver HFODD [43, [43] ca- 
pable of treating simultaneously all possible collective de- 
grees of freedom that might appear on the way to fission. 
In the present work, we adopted the HFODD solver to the 
FT-HFB and FT-HF-^BCS frameworks along the lines of 
SecIH 



IV. RESULTS AND ANALYSIS 

The main objective of this study is to provide a micro- 
scopic description of fission of excited nuclei, based on 
the nuclear DFT. To this end, we solve the constrained 
FT-HF+BCS problem along a collective path defined by 
a mass quadrupole moment Q2o- At each value of Q2o^ 
self-consistent equations are solved, whereupon the total 
energy of the system is always minimized with respect 
to all remaining shape parameters. Along the optimum 
path found in this way, axial and mirror symmetries can 
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Q20 (b) 

FIG. 1: (Color online) Fission pathways of ^''"Pu (top) and 
^^^Fm (bottom) as functions of the mass quadrupole moment 
Q20 at different values of the ground-state temperature kTg.s. 
(marked by numbers, in MeV). Along the minimum-energy 
pathways ("min", solid lines), all self-consistent mean- field 
symmetries can be broken. To illustrate the corresponding 
energy gain, the axial, reflection-symmetric energy curves are 
also shown ("sym", dashed lines). The energy curves have 
been normalized to zero at the ground-state minimum. The 
values of kTg.s.= 1, 1.5, and 2 MeV correspond to excitation 
energies of 13.82, 36.79, and 70.88 MeV for ^'"'Pu, and 14.93, 
39.20, and 75.16 MeV (not shown) for ^'^^Fm. 



be broken, i.e., the multipole moments Q22 and/or Q30 
may be nonzero. Figure [1] shows the fission pathways 
for 240pu and ^ssFm. The former nucleus is known to 
fission asymmetrically while the later one is on the edge 
of the transition from asymmetric to symmetric fission 
EHl- It is, therefore, expected that the fission path- 
ways of these two nuclei would evolve somewhat differ- 
ently with increasing excitation energy. 

For ^■*°Pu, the optimal fission pathway at zero tem- 
perature exhibits the familiar two-humped structure. 
At kTg,s=l.OMeV {E*^13.82MeV), both saddle points 
are reduced by 2-2.5 MeV. The isentropic barriers are 
rapidly quenched with E*, and they become very small 
at kTg,s =2MeV (£'*=70.88MeV) due to the thermal 
melting of shell effects. In order to assess the impact 
of triaxiality on the first, and mirror asymmetry on the 
second saddle point, we computed the axial reflection- 
symmetric energy curve for ^'^'^Pu (marked as "sym" in 
Fig. [1]). The non-axial (Q22) and reflection asymmet- 
ric {Q30) moments along the optimal fission pathway are 
shown in Fig. [21 The energy gain on the first barrier 
due to triaxiality, quite appreciable at r=0, becomes 
practically negligible at fcTg.s.=1.5MeV while the cor- 



responding quadrupole moment Q22 is nonzero even at 
/cTg.s.=2 MeV. This indicates that at large excitation en- 
ergies the energy surface of ^""^Pu becomes very soft in 
the triaxial direction. 

A similar conclusion can be drawn for the reflection 
asymmetric degree of freedom Q^q and its impact on the 
outer barrier. Experimentally, there is clear evidence for 
a transition from asymmetric to symmetric fission with 
excitation energy [23|. The results displayed in Fig. [T] 
are consistent with the observed change in the pattern 
of fission yields. Indeed, at fcTg.s.=2 MeV the calculated 
optimal fission pathway shows a very weak octupole ef- 
fect. 




Q20 (b) 

FIG. 2: (Color online) Variation of non-axial (Q22, top) and 
reflection asymmetric (Q30, bottom) mass moments as a func- 
tion of Q20 and temperature (indicated in MeV) for ^^"Pu. It 
is seen that triaxiality and reflection asymmetry persist to 
kTg, s. =2MeV . However, as indicated in Fig. [1] their impact 
on the total energy is negligible at the largest temperatures 
considered. 

To further explore the transition from asymmetric to 
symmetric fission, we now consider ^^^Fm. In the heavy 
Fm isotopes, a sharp transition has been observed [4J| 
from an asymmetric mass division of spontaneous fission 
products in ^'"'^Fm to a symmetric mass split in ^^^Fm. 
As seen in Fig. [1] and discussed in detail in Ref. [46| . 
at Tg,s,=0 the second barrier along the symmetric fission 
pathway is very broad as compared to the asymmetric 
case, and this explains the asymmetric distribution of 
fission products observed experimentally. However, at 
kTg,s—l-5MeV, the symmetric pathway becomes close 
in energy to the asymmetric one. This indicates that 
competition between asymmetric and symmetric fission 
is expected to occur in ^^^Fm at lower excitation energies 
than in ^^°Pu. 

Wc would now like to address the important question 
of the synthesis of superheavy elements in heavy-ion fu- 
sion reactions. It has already been mentioned that the 
crucial quanity in the synthesis is the survival probablity, 
which depends on the quenching of the fission barrier 
height with E* . In order to obtain a better understand- 
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ing of how the shell effects impact the E* dependence 
of the first saddle point of superheavy nuclei, we per- 
formed systematic FT-HFB calculations for 48 even-even 
nuclei with 110<Z<124 and 166<iV<188. A sample re- 
sult illustrating our methodology is displayed in Fig. [3] 
for Z=112, 118, and 124. 




20 40 60 80 

Excitation Energy (MeV) 

FIG. 3: (Color online) Predicted excitation energy depen- 
dence of barrier heights of even-even superheavy elements 
with Z=112, 118, and 124. 

The dependence of a fission barrier {Eb) on E* is 
usually approximated by a phcnomcnological expression 

Mm 



Eb oc e-''"^', 



(7) 



where the barrier damping parameter 70 characterizes 
the rate of the barrier quenching with excitation energy. 
It is clearly seen from Fig. [3] that the ansatz ([7]) well de- 
scribes the FT-HF-FBCS results and the parameter 
can be meaningfully extracted for every nucleus. This 
is in spite of the fact that many physical effects impact 
Eb-vs-E* dependence. (In addition to a direct depen- 
dence of Eb on entropy, significant contributions come 
from self-consistent variations of nuclear mean fields with 
S, most notably the gradual decrease of triaxiality. The 
quenching of the pairing energy does not impact the ex- 
tracted values of jd as the \ow-E* part of Eb was not 
considered when extracting the slope oi In Eb-) When in- 
specting Fig. [21 one can notice rather dramatic isotonic 
variations of the damping rate for Z—\12. As discussed 
in Ref. [23|, in the isentropic picture, the observed pat- 
tern can be attributed to the higher temperature of the 
lowest minimum as compared to that of saddle point. 

The survey of 7^^ obtained in this work, shown in 
Fig m nicely illustrates the appreciable particle number 
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FIG. 4: (Color online) Inverse barrier damping parameter 7^^ 
extracted from our FT-HF+BCS calculations for 48 even-even 
superheavy nuclei with 110<.Z<124 and 166<Af<188. 



dependence of barrier damping. The maximum of 7^^ is 
predicted for iV=176 and 178, while for iV=166 and 168 
7^^ is fairly small, indicating a rapid decrease of bar- 
rier heights with E* around ^^''112, i.e., in the region of 
deformed superheavy nuclei stabilized by the deformed 
subshell closure N=1Q2 [1, For heavier systems with 
Z=122 and 124, the largest barrier damping effect is ex- 
pected around 7V=182 and 184, i.e., in the region of the 
enhanced shell stability around the expected spherical 
A^=184 magic gap [1, @, 0]- The strong dependence of 
the barrier damping parameter on N and Z indicates the 
importance of shell effects when modeling the formation 
of superheavy elements. 



V. SUMMARY 

In conclusion, we performed systematic self-consistent 
calculations of thermal fission barriers of superheavy nu- 
clei based on the FT-HF+BCS extension of the solver 
HFODD that is capable of describing arbitrary shapes 
free from self-consistent symmetry constraints. Our sur- 
vey of the fission barrier damping parameter demon- 
strates the existence of strong shell effects on 7d. In 
particular, the fastest decrease of fission barriers with ex- 
citation energy is predicted for deformed nuclei around 
A^=164 and spherical nuclei around A^=184 that are 
strongly stabilized by g.s. shell effects. On the other 
hand, for the transitional nuclei around Af=176, the bar- 
rier damping is relatively weak. The particle-number de- 
pendence of 7£i shown in Fig. |4] is expected to impact 
the survival probability of the superheavy compound nu- 
clei produced in heavy-ion fusion experiments; we hope 
that the values of the damping parameter obtained here 
can be useful in guiding future theoretical work on the 
production of superheavy nuclei. 
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We also studied the quenching of triaxial and reflection 
asymmetric deformations with excitation energy. For nu- 
clei ^^°Pu and ^^^Fm, which exhibit asymmetric sponta- 
neous fission, the FT-HF+BCS theory predicts a tran- 
sition to symmetric fission at higher excitation energies. 
Finally, the thermal quenching of triaxiality at the first 
saddle point provides a significant contribution to 70. 
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